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We find the different spatial chaos in a one-dimensional attractive Bose-Einstein condensate inter- 
acting with a Gaussian-like laser barrier and perturbed by a weak optical lattice. For the low laser 
barrier the chaotic regions of parameters are demonstrated and the chaotic and regular states are 
illustrated numerically. In the high barrier case, the bounded perturbed solutions which describe 
a set of discrete chaotic states are constructed for the discrete barrier heights and magic numbers 
of condensed atoms. The chaotic density profiles are exhibited numerically for the lowest quantum 
number, and the analytically bounded but numerically unbounded Gaussian-like configurations are 
confirmed. It is shown that the chaotic wave packets can be controlled experimentally by adjusting 
the laser barrier potential. 
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I. INTRODUCTION 

It is well-known that chaos in a nonlinear system not 
only may play a destructive role^but also has many prac- 
tical and dramatic applications [lj . Chaos has been thor- 
oughly studied during the last century in many differ- 
ent fields of physics. Very recently it has been recog- 
nized that existence of chaos is also possible in the Bose- 
Einstein condensates (BECs) described by the Gross- 
Pitaevskii equation (GPE) picture [2| and in the dis- 
cretized systems describingtrapped BECs within the 
space-mode approximation [J]. The temporal chaos was 
revealed in the time evolutions of BECs trapped in a 
double- well potential The spatial chaos with spa- 
tially disordered configurations was investigated for the 
stationary states of the BECs held in an optical lattice 
5]. The spatiotemporal chaos in BECs interacting with 
different potentials has also been found @. 

The mean-field stationary states of a BEC are dom- 
inated by the time-independent one-dimensional (ID) 
GPE M,]!]. For a BEC in ground state without cur- 
rent [3] the GPE is a real equation and can be identi- 
cal with the celebrated Duffing equation 0, [l(| by using 
time instead of spatial coordinate. Particularly, when 
such a GPE is perturbed by a weak periodic potential, 
the Smale-horseshoe chaos may appear for a certain pa- 
rameter region of the extended dynamical system [f| [ic| . 
The Melnikov chaos criterion gives the chaotic parame- 
ter region in which the perturbation parameters are al- 
lowed to vary their values continuously [III, [H, 0, [3] • 
The Gaussian-like barrier potentials can be realized by a 
sharply focused laser beam in the experiment [TH ] , which 
have been applied to investigate the shock-wave forma- 
tion in BECs the nonlinear resonant transport (l7j 
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and deterministic chaos [l8| of BECs. Recently, using 
external fields to control quantum states of BECs has 
become an important physical motivation [19j . 

When a BEC is created initially in a time-independent 
optical lattice, the stationary states of the GPE are de- 
termined by the boundary conditions and are adjusted 
by the system parameters. The different boundary condi- 
tions may be established in a practical experiment, which 
cannot be set accurately. In the non-chaotic regimes, a 
small change of the boundary conditions and/or system 
parameters brings the BEC state only a small correc- 
tion which can be neglected in a good approximation. 
In the chaotic regimes, however, the stationary state de- 
pends on the conditions and parameters sensitively. The 
sensitivity means that a small change of the conditions 
and/or parameters may cause a great difference which is 
not negligible. For example, the periodic configuration 
of BEC density is changed to the aperiodic and irregular 
one. It is important for the application purpose to predict 
the bounded states and to manipulate the corresponding 
density distributions which govern the beam profile of an 
atom laser extracted from the BEC [2fJ. Therefore, in- 
vestigating the spatial chaos and its control is necessary 
and interesting for the considered BEC system. 

The main aim of this paper is to present an analytical 
evidence of a different type of spatial chaos which can 
be defined as the discrete chaotic states, and to estab- 
lish a method for controlling the chaotic states. By the 
discrete states we mean a denumerable set of bounded so- 
lutions in which any solution is one-to-one with a value 
in a discrete set of the parameter values. If the discrete 
states meet the Melnikov chaos criterion, we call them the 
discrete chaotic states. By using a laser beam modeled 
by the tanh-squared-shaped barrier potential [2l| which 
is known as the Rosen- Morse potential [13], we demon- 
strate the existence of spatial chaos in the BEC held in 
a weak optical lattice. The chaotic regions of parameters 
are exhibited and the regular and disordered configura- 
tions of the BEC are illustrated. It is shown that the 
width and site of the strong barrier potential confine the 
width and site of the BEC wave-packet, and a denumer- 
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able set of the barrier height values corresponds to the 
discrete chaotic states and magic numbers of condensed 
atoms. Thus the possible chaotic states can be controlled 
by adjusting the width, site and height of the laser barrier 
experimentally. 
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II. CHAOTIC AND REGULAR STATES FOR 
THE LOW LASER BARRIER 

For the considered BEC system with transverse wave 
function being in ground state of a harmonic oscillator of 
frequency w r , the governing time-independent quasi-lD 
GPE reads 



2m 



(1) 



where m is the atomic mass, \i is the chemical potential, 
and g' ld = gomuj r / (2irh) = 2ftw r a s denotes the quasi-lD 
atom-atom interaction intensity with a s being the s-wave 
scattering length. Hereafter, by ip xx we mean the second 
derivative of ij) with respect to x. The external potential 
V'{x) — — Vo tanh 2 [/3(x — x c )] + V\ sin 2 kx contains the 
longitudinal barrier potential of strength Vq > 0, width 
and center site x c , and the perturbed lattice poten- 
tial with Vi and k being the intensity and wave vector. 
The former as a Gaussian-like potential can be formed by 
a sharply focused laser beam in the experiment [l5[ , and 
the latter is a laser standing wave. Taking and f3 as 
the units of coordinate x and density \ip\ 2 , and normaliz- 
ing the potential strengths Vq , V\ and chemical potential 
/i by using Ep = h 2 f3 2 /m, Eq. (1) becomes the dimen- 
sionless equation 



-xP xx + [V{x)+g ld \^\ 2 ]^^ ^. 



(2) 



Here the interaction intensity is reduced to 
2hw r a s (3 / Ep — 2a s /(/3a 2 ) with a 



9id = 

y/H/(muj r ) being 



the transverse harmonic oscillator length, and the poten- 
tial gets the form 



V(x) = -Vq tanh 2 (x - x c ) + Vi sin 2 kx 



(3) 



with k measured in /3. 

We are interested in the real solution of GPE (2), 
which makes the GPE the perimetrically perturbed Duff- 
ing equation (Toj in the spatial evolution and for the weak 
potential. It is well known that existence of the periodic 
perturbation is necessary for the appearance of chaos in 
the Duffing system [H El El- When negative inter- 
action and negative chemical potential are taken, in the 
absence of external potential the system has the well- 
known homoclinic (separatrix) solution [T3, EH, EH) Ell 
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FIG. 1: (a) The constant F as a function of Co — x c for pa- 
rameters fi = —2 and Vo = 0.2. (b) The boundaries of the 
chaotic regions for fi — —2, and F = 0.005 (dashed curve), 
F — 0.01 (solid curve) and F = 0.015 (doted curve). 



where Co is an arbitrary constant adjusted by the bound- 
ary conditions at the boundary x = xo- For the BEC 
system governed by Eq. (2) the constant Co cannot be 
determined experimentally, because of the undetectable 
■00 (^o ) ■ The presence of the weak external potential leads 
to the Melnokov function fiR EH, E3, El 



M(co) 



gid 



F 



24> 0x {x)V{x)ilJo(x)dx 
fc 2 7rVi sin(2fcc ) 



2/i sinh(fc7r/ \J— 2/1) 



(5) 



for < Vo <SC 1 and |Vi| -C 1, where ipox denotes the first 
derivative of ipo with respect to x, constant F from the 
barrier potential reads 



4V e 



2A 



lMl(e 2A -l)^ 



[(3 + 2A + 8Ae + (2A - 3)e 4A )](6) 



with A — Co — x c . The Melnokov function measures the 
distance between the stable and unstable manifolds in the 
Poincare section of the equivalent phase space {tjj, ipx)- 
For some Co values if the Melnikov function has a simple 
zero, the locally stable and unstable manifolds intersect 
transversally such that the Smale-horseshoe chaos exists 
in the Poincare map [l(| EHj E2, EH- The possibility 
of M{co) = results in the chaotic region of parameter 
space 
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When parameters are taken in the chaotic region, the 
Melnikov function has zero points and the stable and un- 
stable manifolds in the Poincare section may intersect 
that leads to the Smale-horseshoe chaos. It is possible 
that the regular orbits exist for both the chaotic and 
non-chaotic regions. The chaotic and regular orbits in 
the chaotic region depend on the different boundary con- 
ditions respectively. 

As can be seen from Eq. (7), for any negative chemical 
potential /i < and any barrier potential strength in the 
region < Vb <C 1, the chaotic region depends on con- 
stant F in the plane of parameters V\ versus k. The F 
is determined by the parameters Vo, [i and Co — x c with 
the potential strength Vq and site x c being adjustable. 
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FIG. 2: The Poincare section on the equivalent phase space 
(if), ipx) from Eq. (2) for the given parameters and boundary 
condition. 



In Fig. la we show F as a function of cq — x c for /i = —2 
and V = 0.2 by using the MATHEMATICA code. From 
this figure it can be observed that has a maximum 
\F\ = 0.015 and a minimum \F\ = 0. The former cor- 
responds to the minimal chaotic region of Eq. (7), and 
the latter is associated with the maximal chaotic region 
|Vi | > 0. Taking fi = -2 and F = 0.005, 0.01, 0.015 as- 
sociated with three different Cq values respectively, from 
Eq. (7) we plot the boundary curves of the chaotic re- 
gions as the dashed curve, solid curve and doted curve 
of Fig. lb. The corresponding chaotic regions are above 
these curves respectively. The minimal chaotic region 
above the curve of |F| = 0.015 is certainly chaotic re- 
gion for arbitrary Co value. But the other chaotic regions 
are related to the corresponding boundary conditions, 
through the constant F(cq). 

A useful way of analyzing chaotic motion is to look at 
what is called the Poincare section, which is a discrete set 
of the phase space points at every period of the periodic 
potential, i.e. at x = 2ir/k, 4ir/k, 6ir/k, ■ ■ ■ . Taking the 
parameters fx = —2, g ld = —1, V Q = 0.2, V\ — 0.2, k — 
1.5, x c — 1 and the approximation [\l>(xo),^ x (x(yj\ — 
[V>(1000Q),Vx(10000)] = (0.00001,0.00001) to the exper- 
imentally possible boundary condition [^(oo), ^(oo)] = 
(0,0), from Eq. (2) we numerically plot the Poincare 
section on the equivalent phase space (jp, i/) x ) and find 
the chaotic trajectory as in Fig. 2. Here the lattice 
strength V% and wave vector k are evaluated in the min- 
imal chaotic region of Fig. lb. For the same parame- 
ters of Fig. 2 from Eqs. (2) and (3) the potential and 
chaotic state functions are plotted as in Figs. 3a and 
3b respectively. From Fig. 3a we can see the profile 
of the combined potential between the barrier potential 
and periodic lattice. In Fig. 3b we exhibit the aperi- 
odicity and irregularity of the chaotic macroscopic wave 
function corresponding to Fig. 2 numerically. In order 
to confirm the sensitive dependence of chaotic system on 
the boundary conditions, we change only the boundary 
condition as [^(10000), ^ x (10000)] = (0, 0.00001) to plot 
the wave function. This small change leads the irregular 
curve in Fig. 3b to the periodic one in Fig. 3c. When 
the lattice strength is decreased to V\ = 0.005 and the 
other parameters are kept, from Fig. lb we observe that 
the parameter value is outside the given chaotic region. 
After changing V\ from 0.2 to 0.005, Figs. 3a, 3b and 



3c are changed to Figs. 4a, 4b and 4c respectively. Fig- 
ure 4a displays the weak periodic potential compared to 
the laser barrier. In Fig. 4b and 4c we illustrate that 
in the considered parameter region the wave functions 
are periodic for the given boundary conditions. It is in- 
teresting noting that the regular wave functions in Figs. 
4b have two different periods and two different ampli- 
tudes in both sides of the laser barrier. This means that 
the atomic number J s \%j}(x)\ 2 dx is different for the in- 
tegration region £ of different side. The periodicity is 
varied with the change of the boundary conditions from 
h/>(10000), V'z(lOOOO)] = (0.00001,0.00001) of Fig. 4b to 
[^(10000), ^(10000)] = (0, 0.00001) of Fig. 4c. Differ- 
ing from Fig. 4b, in Fig. 4c the period of wave function in 
both sides of barrier is the sameness and the smaller one 
of amplitudes is enlarged compared to that of Fig. 4b. 
The results display the different profiles of macroscopic 
quantum states and reveal that the existence of chaos 
means the sensitive dependence of the BEC system on 
the boundary conditions and parameters. 



III. DISCRETE CHAOTIC STATES FOR THE 
HIGH LASER BARRIER 

The chaotic region of Eq. (7) is based on the perturba- 
tion theory [TTL [l2j so that it is valid only for very small 
potential strengths Vo and V\. When the strength V) of 
the barrier potential is continuously increased to strong 
enough, e.g. Vq > 1, it can no longer be treated as a part 
of perturbations. In this case we require to reconsider the 
perturbation problem of the stationary states. Applying 
the well-known Rayleigh-Schrodinger expansions [231 ] 

ip = i/} + ipi, M = Mo + Mi for |t/>i|, |Vi|<l(8) 

to Eq. (2) of real ip, we have the leading order and the 
first order equations as 

--ipoxx - [Vo tanh 2 (x - x c ) - SidV'olV'o = Mo^o,( 9 ) 



■■^ipixx - [V tanh 2 (x - x c ) - 3g ld tp^ + ^ ]V>i 
= (fj,i — Vi sin 2 kx)ijjo(x). 



(10) 



Noticing that Eq. (9) has many special solutions for the 
fixed values of Vo, gu, x c and different /iq values. Only 
the homoclinic solution is related to the Melnikov's chaos 
and the other solutions are associated with the regular 
states of Eq. (2). Here we are interested in the chaos and 
only consider the homoclinic solution thereby. It can be 
easily proved that the homoclinic solution of Eq. (9) has 
the form 



ipa = \ ° + sech(a; - x c ) for ^ = -V - -■ (11) 

Differing from Eq. (4), Eq. (11) describes a wave packet 
whose height and width are adjusted by the potential 
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FIG. 3: (a) The potential function of Eq. (3) and (b) the aperiodic chaotic state of Eq. (2) for the same parameters and 
boundary condition with Fig. 2. (c) When the value of ?/>(10000) is changed from 0.00001 to 0, we get the periodic wave 
function. 




FIG. 4: The correspondences of Fig. 3 after the change of parameter Vi from 0.2 to 0.005. 



intensity Vq and width implied in the unit of x. 
Substituting Eq. (11) into Eq. (10) yields the non- 
homogeneous equation 



q 1 
(2V" + 3)secli (ir-Xc) - - 



-f(x) = {fix - Vi sin 2 kx)ip (x). 



(12) 



The corresponding homogeneous equation for / = is 
a well-known Schrodinger one with trapping potential 
— (2Vq + 3)sech 2 (a; — x c ) and eigenenergy E = —1/2. 
Given two linearly independent solutions of the homo- 
geneous equation as and i)>'{ = ip[ J (ip[)~ 2 dx, the ex- 
act general solution of non-homogeneous Eq. (12) can be 
written in the integral form [24] 

f-X f-X 

i> x = 2^4' / ^'J(x)dx - 2^i / i>'{f(x)dx, (13) 

J A JB 

where A and B are arbitrary constants determined by 
the boundary and normalization conditions. This so- 
lution can be directly proved by comparing the second 
derivative tp\ xx from Eq. (13) with that in Eq. (12). 

Boundedness of the perturbed correction ipi is the 
physical requirement, which depends on the bounded ip^. 
In order to seek such a ip[, we set [23| 

ip[ = [sech(a; — x c )] 2X u(z), z — — sinh 2 (a; — x c ), 
A = [v/8(2Vb + 3) + l-l]/4. (14) 

Inserting Eq. (14) into the homogeneous part of Eq. (12) 
with / = produces the hypergeometric equation 

z(l - z)u zz + [0.5 - (a + b + l)z]u z - abu = 0, (15) 

where a = 0.5 — A, b — —0.5 — A. Its two linear indepen- 
dent solutions with finite terms read [23j 

< = ^(0.5 — A, -0.5 - A, 0.5, z) for A = 0.5 + n; 

u° n = y/\z\F(l — A, —A, 1.5, z) for A = 1 + n. (16) 



Here F(a, b, c, z) is the hypergeometric function, u e n and 
u° n with n = 0, 1, 2, • • • denote even and odd functions of 
(x — x c ) respectively. Combining Eq. (16) with Eq. (14) 
we arrive at the bounded solutions 

= [sech(x - x c )] 1+2n u^ V = [(3 + 4n) 2 - 25]/16, 
ip[° n = [sech{x - x c )] 2+2n u° n , V Q = [(5 + An) 2 - 25]/16 



for Vq > 0, ra = l,2,- 



(17) 



Noting that tp[ of Eq. (17) tends to zero, and tp" = 
ip[ J (i/)^) -2 ^ of Eq. (13) is infinity at x = ±oo. Thus 
the second term of Eq. (13) is in the form of zero multi- 
plying infinity at x — ±oo, so we can use the l'Hospital 
rule to calculate the limit and to prove the boundedness 
of this term. For the first term of Eq. (13) we have to 
establish the boundedness condition 



lim 

x— *-±oo 



V>i(Fisin 2 fcz;- u 1 )ip o dx = 0. (18) 



The necessity of Eq. (18) is obvious for the bound- 
edness of Eq. (13), because of the unboundedness of 
tp". Under condition (18) we can apply the PHopital 
rule to the both terms of Eq. (13), obtaining [24| 
linij^ioo ipi = 2]im x —>± 00 f(x) — 0. This limit implies 
that Eq. (18) is also sufficient and the obtained macro- 
scopic wave function satisfies the usual boundary condi- 
tion ijj(±oo) = ^o(±oo) + -0i(±oo) = 0. Noticing the 

correspondence between f(x) and e^ 1 (x) in Eq. (9) of 
the second article of Ref. [2J|, the above proof of suffi- 
ciency is clear. 

The integration of the first term in Eq. (13) is in- 
solvable and cannot be expressed by finite elementary 
functions. Hence, in the numerical computation based 
on Eq. (13), small deviation from the exact value of the 
integration satisfying condition (18) is avoidable. The 
small deviation will be amplified exponentially fast by 
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the unbounded function if>i(x) until infinity as x — > ±00 
that exhibits the numerical instability. The analytical 
insolvability and numerical instability can cause the un- 
predictable chaotic behavior [l4[ • The difference I + — I_ 
of the integration in Eq. (18) is similar to the Melnikov 
function of Eq. (5), hence /+—/_ = can be called the 
generalized Melnikov criterion for chaos. In fact, given 
Eq. (18), the undetermined form -0"(±oo) x I± = 00 x 
appears in Eq. (13) as |a;| tending to infinity, which 
leads to a new feature as the analytical boundedness but 
numerical unboundcdncss, namely the evidenced incom- 
putability and unpredictability of the chaotic behavior 
pL4| . Therefore, under the condition (18) the solution 
ip(x) = ipo(%) + i>i(x) in terms of Eqs. (11) and (13) 
is called the chaotic solution • If the zero boundary 
condition [^(±00), ^(ioo)] = (0,0) is required theoret- 
ically, the uniqueness theorem infers the chaotic solution 
to be the unique one of the system. On the other hand, 
from the formula of the energy functional @, Q , 



H = 



- 7V 



1 + V(f) 4 



ipd 3 x 



^9id\^\ 4 



d 3 :i 



(19) 



we know that unlike the unbounded solution with 
(^+(±00)1 = 00, the analytically bounded solution with 
-0 + (±oo) — is associated with the finite energy func- 
tional and may be metastable thereby Q. Although the 
chaotic solution is not very stable, due to the sensitive 
dependence on the parameters and boundary conditions, 
it may also be metastable compared to the analytically 
unbounded solution. Particularly, these bounded solu- 
tions are valid only for the discrete Vq values of Eq. 
(17). This means the corresponding analytically bounded 
chaotic states to be discrete with the increase of the bar- 
rier height. 

The above-mentioned results imply that when the bar- 
rier potential is strong enough, its strength values must 
be discrete for the bounded perturbed solutions. For the 
discrete Vq = Vbn values the leading number-density ipQ 
is proportional to Vbn and the leading chemical potential 
is given as ^ n = — \ — Vo n by Eq. (11), the both are 
also discrete. The parameters V\,k,x c and g\d can vary 
their values continually in a certain parameter regions. 
Given a set of values of V±,k,x c , the first correction \i\ 
is determined by the boundedness condition of Eq. (18). 
In Eq. (10) the discrete chemical potential /i ss /j,Q n + /ii 
is equivalent to the energy of a Schrodinger system. In 
quantum mechanics [HI], it is known that the bounded- 
ness of wave function may lead the energy to take discrete 
values. Mathematically, the relationship between the dis- 
crete values of potential strength Vo = Vbn and the ex- 
actly bounded solutions of Eq. (12) agrees qualitatively 
with that of a 2D Coulomb correlated system (2(| , where 
the Schrodinger equation is exactly solvable only for a de- 
numerably infinite set of values of magnetic strength (or 



the corresponding oscillator frequency). Physically, we 
well know that for a 2D electron gas in a semiconductor 
heterojunction the integral and fractional quantum Hall 
plateaus are associated with the discrete set of values of 
magnetic strength [27 1. 

We now investigate the physical effect of the discrete 
laser strength Vo — Vbn on the considered BEC system. 
Applying Eq. (11) to the normalization condition yields 
the number of condensed atoms N n w J \ip 0n \ 2 dx = 2(1 + 
K)n)/|<7id| = (1 + Vo n )(3af,/\a s \ for the metastable states 
given by Eq. (8) with Eqs. (11) and (17), that results in 
the relation 



N n \a s \ « {1 + V 0n )pa 2 r 



(20) 



with Vbn given in Eq. (17). Here the special value N n can 
be called the magic numbers of the macroscopic many- 
body system keeping in the metastable states. Differing 
from the magic numbers of the microscopic many-body 
system (e.g. atomic nucleus), N n denotes some approx- 
imate values, because of the approximation N ± 1 w N 
in the mean-field theory of macroscopic many-body sys- 
tem 0, H| . For a harmonically confined BEC system, the 
supercritical number N cr of condensed atoms obeys Q 
A^cr|&s| = 0.575a/! O with a^ being the 3D harmonic os- 
cillator length. The magic number N n may exceed the 
supercritical number N cr by increasing the laser strength 
Vbn and/or decreasing the laser barrier width The 
approximate magic numbers of the considered many- 
body system warrants experimental investigation. 

Let us take the simplest even solution of [x — x c ) with 
quantum number n = 1 as an example to show the fea- 
ture of the chaotic solutions. From Eqs. (17) and (16) 
such a solution is derived as 



ip[ n = tpxi = sech 3 y(l - 4sinh 2 y) 



(21) 



for y — x — x c and Vbi = 3/2, /101 = —1/2 — Vbi = —2. 
Obviously, this solution tends to zero as x — > ±00. The 
corresponding unbounded solution reads 



1 

= — sech y (36y — 24y cosh 2 y + 28 sinh 2y — sinh iy) 
64 

(22) 

in which the term sech 3 y sinh 4y tends to ±00 and the 
other terms tend to zero as x — > ±00. Applying Eqs. 
(21) and (22) to Eq. (13), the exact general solution of 
Eq. (12) becomes 



ru = 2iPn / V>i e i/(^ - 2^1 / ^f(x)dx, (23) 

J A J B 

where f{x) = — (^1 — V\ sin 2 kx)ipo{x) is equal to zero 
at x = ±00, because of ^(±00) = 0. Applying the 
l'Hopital rule to Eq. (23), we easily verify its bounded- 
ness [24j . through the limit \mi x ^± 00 tpl^x) = for the 
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fii obeying Eq. (18) accurately. Inserting ip£ and Eq. 
(11) into Eq. (18), from the generalized Melnikov chaos 
criterion I + — J_ = one derives 

Mi = 0.5ttVi cos(2kx c )k(5k 2 - l)csch(/c7r) (24) 

which can be adjusted by the laser site x c and has a max- 
imum and a minimum at cos(2fca; c ) = ±1 respectively. 

In order to obtain the bounded numerical solution of 
Eq. (23), the parameter /ii must obey Eq. (24). How- 
ever, in any numerical computation, for a set of fixed 
parameters V\, k, x c it is impossible to take the value 
of fii accurately, because of the irrational 7r with infi- 
nite sequence of digits in Eq. (24). This implies small 
deviation from the accurate boundedness condition (18) 
and the small deviation will lead the numerical solution 
of Eq. (23) to be exponentially amplified by the un- 
bounded function until infinity as x — ► ±oo. So the 
analytically bounded chaotic solution (23) is numerically 
unbounded and uncomputable for sufficiently large |x| 
values [l4j. For a small \x\ value -011 is finite and Eq. 
(23) is certainly bounded. At x = ±oo the boundedness 
condition (18) and PHopital rule lead Eq. (23) to zero an- 
alytically. The unpredictability of chaotic solution (23) 
may occur only near the spatial range \y\ — |x — a; c | G 
(|y s |, oo), where y s — x s — x c can be estimated through 
the starting point of the numerical incomputability after 
which the solution tends to infinity rapidly. In such a 
spatial range, the chaotic region of atomic density may 
be |V(z/)| 2 6 (2|^b(±oo)Vfi(±oo)|, 2|^,(i/.)Vfi(tf.)|) = 
(0, 2\iJjo(y s )ipf 1 (y s )\ with width S(y) tending to zero as 
the increase of \y\ value. The maximal width reads 
6(y s ) w 2\ipo(y s )ipf 1 (y s ) \ which is in order of perturbation 
Vi, since \ip{y)\ 2 equates \ip (y) + ^n(y)| 2 ~ \My)\ 2 + 
2|V'o(2/)'0ii(j/)l aim l^oiy)] 2 is predictable for any y. The 
effective first-order correction to the Gaussian-like pro- 
file is analytical bounded, which can be obtained by cut- 
ting the infinity from the numerical solution of Eq. (23). 
These will be illustrated numerically as follows. 

As an instance, setting the parameters Vq = Vol = 
3/2, Vi = 0.05, k = 1.5, g u = -l,/x = -2 and the 
boundary condition which is equivalent to A = — oo, B = 
0, from Eqs. (11) and (23) we plot the chaotic atomic 
density \ip\ 2 = (?Ao + V'n) 2 as in Fig. 5a. Here the solid 
and dashed curves correspond to x c = l,fii = —0.02148 
and x c = 2,fii = 0.02083 respectively, which satisfy the 
generalized Melnikov chaos criterion (18) and (24) ap- 
proximately. The dashed curve has approximate shape 
with the solid one and can be regarded as the latter after 
a translation of distance 1. The numerically unbounded 
first corrections are uncomputable for sufficiently large 
\y\ = \x — x c \ values and the starting points of the incom- 
putability are shown to be about y = ±y s « ±2 after 
which the atomic densities may be irregular and tend to 
infinity rapidly. By using the wide-black curves instead 
of the infinity in range \x — x c \ > 2 of Fig. 5a, we ob- 
tain the Gaussian-like wave packets as in Fig. 5b which 
describe the analytically bounded atomic density better. 
The wide-black parts are the sketch maps of the chaotic 





FIG. 5: (a) The chaotic density profiles of atomic number 
for x c = 1 (solid curve) and x c = 2 (dashed curve), (b) 
The analytically bounded density profiles from Fig. 5a by 
replacing the parts of \x — x c \ > 2 with the sketch maps of 
the chaotic density regions. 



regions of density distributions, whose width varies from 
the maximal value S(y s ) ~ V\ to minimal one <5(±oo) = 0. 
In the chaotic regions of density, the atomic density is un- 
predictable. The effective first corrections in the range 
x G (—0.5, 3) are exhibited by the inset of Fig. 5b, which 
are plotted from Eq. (23) for the range \y\ < 2 and 
the parameters adopted in Fig. 5a. It should be em- 
phasized that the analytically bounded chaotic states are 
discrete and can be manipulated experimentally by tak- 
ing the barrier heights Vo« in Eq. (17) discontinuously 
and adjusting the barrier site x c continuously. Partic- 
ularly, by increasing x c adiabatically [28j, we can move 
the Gaussian-like wave packets slowly for the purpose of 
BEC transport 



IV. CONCLUSIONS AND DISCUSSIONS 

We have investigated the spatial structure of the ID 
attractive BEC interacting with a tanh-squared-shaped 
laser barrier potential and perturbed by a weak laser 
standing wave. The existence of the Smale-horseshoe 
chaos is demonstrated and the Melnikov chaotic regions 
of parameter space are displayed. In the low laser bar- 
rier case, the aperiodic chaotic states and periodic regu- 
lar states are illustrated numerically. For the sufficiently 
strong barrier potential a set of discrete chaotic solu- 
tions is constructed formally. Any chaotic solution is 
the combination of a Gaussian-like wave-packet with the 
corresponding perturbed correction. The discrete chaotic 
solutions are analytically bounded only for the discrete 
barrier height values and special magic numbers of con- 
densed atoms. The density profiles of BEC in the dis- 
crete chaotic states are investigated numerically for the 
lowest quantum number, and the numerical instability is 
revealed. The Gaussian-like wave could be translated by 
varying the laser-barrier site adiabatically, which is sim- 
ilar to the bright soliton of an attractive BEC with the 
parabolic barrier potential [25| . The periodic structures 
of BEC can be detected by the Bragg scattering of an 
optical probe beam [2^] and the used Gaussian-like po- 
tential can be generated by a sharply focused laser beam 
in the experiments [l5| . Thus the irregular chaotic states 
could be observed and controlled readily with current ex- 
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pcrimental capability. 

The existence of chaos means the sensitive dependence 
of the BEC system on the boundary conditions and pa- 
rameters in chaotic region. The sensitivity causes the 
unpredictability of the spatial distributions of the BEC 
atoms, since the boundary conditions cannot be set ac- 
curately in a real experiment. The above results reveal 
the possible bounded states associated with the spatial 
distributions, and suggest a method to control the irreg- 
ular chaotic states by adjusting the lattice strength and 
laser barrier parameters. 

It is worth noting that the discrete chaotic states may 



appear in many different physical systems with different 
Gaussian-like potentials and may also exist in the tempo- 
ral and spatiotemporal evolutions of the time-dependent 
systems. 
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